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Abstract 

We calculate the x dependence of nuclear shadowing at moderate values of by 
using HERA diffractive data. We show that no decrease of shadowing occurs down to 
very low x (x ~ 10~^). 
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The possibility of studying nuclear shadowing {i.e. the depletion of bound nucleon 
structure functions (-^2^^) with respect to the free nucleon ones (F^)) at HERA ^ has 
prompted new interest on this subject. Different models and interpretations have been 
suggested in the past years to explain nuclear shadowing . A nuclear HERA program 
would allow extending the experimental investigation to very low x, thus offering the 
chance of a deeper understanding of the phenomenon. 

One recent interesting prediction 0] is that, at low x and for moderately large 
{Q^ ~ 10 GeV^) shadowing stops its rise, starts decreasing and eventually vanishes. 
This effect arises from the different low-x behavior of the inclusive and the diffractive 
structure functions. In fact, at x ^ 1, whereas F2{x,Q'^) ~ {l/x)^^^'^\ with a 'hard' 
intercept A(Q^) ^ 0.3 — 0.4, the shadowing term AF^*^ = — F^ , which is propor- 
tional to the diffractive structure function, is dominated by soft physics and has an 
asymptotic x dependence of the form AF|^ ~ (l/x)^^*-^^-*, with A(yU^) ^ 0.1, where /i^ 
is the typical scale of soft interactions. 

While the argument above is rather general and correct, the specific finding of 
Ref. [|], i.e. that at ~ 10 GeV^ shadowing reaches its largest value already at 
X ~ 10^^ and then begins to decrease towards smaller x, is based on some restrictive 
assumptions. We shall call xq the x value where shadowing attains its maximum. The 
value of xq is found in to be almost independent of (for > 5 GeV^) and of 
the atomic mass. The prediction of is that in the region between x = 10~^ and 
X = 10~^ shadowing should be a decreasing function of x towards x — 0. It is evident 
that such a behavior should easily be observable at HERA. 

The results of Ref. rely on two assumptions: 

1) the shadowing term, which is related to the diffractive structure function, is given 
by a simple power-like parameterization 

AFf ~ a;2-2"r{0) 

with aip(O) ~ 1.1; 

2) the mass of the diffractively produced states is kept fixed and equal to Q^. 
Both these assumptions are questionable. The first one is true only asymptotically, 

as X ^ 0, and cannot be used to draw any conclusion about the x-region between 10""^ 
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and 10~^. The second assumption is not valid when the dynamics is dominated by 
the triple pomeron contribution, which provides a nonvanishing high-mass tail to the 
distribution of the diffracted states. 

The purpose of this letter is to carry out a more precise phenomenological analysis 
of nuclear shadowing, by relying on the diffractive structure function measurements of 
HERA. Our main concern will be with the x-dependence of shadowing, rather than 
with its absolute normalization, which is at present unpredictable due to experimental 
and theoretical uncertainties. Our calculation, which does not make use of the assump- 
tions of Ref. 0, leads to the conclusion that the onset of the decrease of shadowing is 
likely to be around or smaller than xq = 10~"^ and thus hardly reachable at HERA. In 
other terms, we predict that no decrease of shadowing will be observed at HERA. 

Let us start from the well-known relation between nuclear shadowing and diffraction 
dissociation established long ago by Gribov ^. In virtual-photon-nucleus scattering 
the nuclear cross section is given by 0, ^ 

^A^r« _ 4. ^ J,,,. J^l^^^ /d^6|*(.,.6)P + . . . , (2) 

where d^a^ / dM'^dt is the 7*A^ diffraction dissociation cross section (M^ being the 
invariant mass of the excited hadronic states), and the longitudinal form factor $, 
which is function of the impact parameter h and of the longitudinal momentum of the 
recoil proton ki = xrriNil + M^/Q^), is related to the nuclear density pa by 

^{kLM= UzpA{b,z)e'^^' . (3) 



The dots in eq. (0) represent the higher rescattering terms, which are non negligeable 
for heavy nuclei. The simplest way to take them approximately into account is ^ to 
introduce an eikonal factor e~°"'=«^*^^^/^ in the integration over the impact parameter in 
(0). T{h) is the nuclear thickness 

/+00 
dzpA{b,z)- (4) 
-oo 

(Teff is the effective cross section for the interaction of the diffracted states with the 
nucleon, given by 



aeff = levr-^ ^ . (5) 
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At small X this cross section turns out to be almost x independent. In practical 
calculations it will be taken as a constant (see below). 

In terms of the inclusive structure functions per nucleon F^'^, and of the diffractive 



structure function 



Q2 dpdQ^dxjpdt 



(6) 



i=0 



where (W is the squared center-of-mass energy of the 7*A^ system) 



_ M' + Q' X 



eq. becomes (k^ = rrij^xp) 

F^{x,Q') = F^{x,Q') 

- 47r^ ldxjpF,''^'\p,Q',xjp,t = 0) J d^be'^^^^^^ Mxjp,b)m 

In principle, the diffractive structure function F^^^^ can be obtained from the exper- 
iment. However, what the present experimental analyses provide is only the diffractive 
structure function integrated over t 

F^'''\(3,Q\xjp) = r^^^d\t\F^''^{(3,Q\xjp,t), (9) 

JO 

with tmax — 0.5 GeV^. 

The ZEUS parametrization for F2^^^ has the form 

= Ax-^^[(3{l-(3) + ^{l-(3f] (10) 

where a = 1.30 ± 0.081^4, A = 0.018 ± 0.001 ± 0.005, C = 0.57 ± 0.12 ± 0.22. The 
exponent a is found to be essentially independent of /3. A more recent preliminary 
analysis [§] gives a smaller value for a: a ~ 1.1 — 1.2 in the range 10 — 20 GeV^. In 
eq. (p!0|) we separated a 'valence' part F^^^^ oc /3(1— /?) and a 'sea' part F^^^^ oc (1— /9)^. 



In the language of the color dipole model [0, the valence corresponds to the lowest 
Fock state (qq) of the virtual photon, whereas the sea corresponds to higher Fock states 
{qqg ■■■)■, which represent the triple pomeron contribution of Regge theory. Notice that 
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factorization breaking effects, predicted in Ref. [11|, which would imply a non-universal 
flux factor for the valence and sea components, are not yet observable. 

In order to derive F^'''^\(3,Q'^,X]p,t = 0) from the measured F^^^\(3,Q'^,X]p), 
eq. (0), we assume a simple peripheral t-dependence of the form 

KTiP,Q^xjp,t) = FZ^^\p,Q^xjp,t = 0)e-(^-=+^"-'°^^)l*l , (11) 

where a'p ^ 0.5 and the slopes B^^s are borrowed from hadron scattering and real 
photoproduction. We use 5^ ~ 6 GeV"^, 5^ ~ 12 GeV~^. 

In eq. (g) F^^^\p,X]p,Q'^) is integrated over xjp between x and 1. However, due to 
the selection of the rapidity gap events, there is an experimental upper cutoff on xjp: 
x']p = 0.04. We checked that in the x region of interest the large-xjp tail neglected in 
the integration is irrelevant. 

For consistency, we used the ZEUS parametrization for F^ix^Q"^) [|l^]. Our pre- 
dictions do not depend on the parametrization adopted for F2 (the same results are 
obtained using MRS(Rl) [0). As for the nuclear part of the calculation, we used a 
Fermi-type nuclear density. 

The two main sources of uncertainty in our calculation are: i) the effective cross 
section (a constant value, specified below, is used); ii) the large experimental error 
on a, the exponent of the so-called pomeron flux, see eq. (0). As we shall see, these 
uncertainties prevent us from predicting the absolute amount of shadowing, although 
they do not affect the qualitative features of the x-dependence of shadowing, that we 
are most interested in. Both the normalization and the x-dependence of shadowing 
depend very little on the other parameters appearing in the calculation. In particular 
we checked that even a large variation of the slopes B^^s changes by no more than few 
percent the predicted shadowing. 

Let us come now to the results. 

In Fig. 1 we plot the ratio F^jF^ for calcium. The diffractive structure function 
used in the calculation is given by the ZEUS parametrization (|10|), which is valid at 
moderately large (Q^ ~ 10 — 30 GeV^). We use the central values of the ZEUS 
fit for the coefficients A and C, quoted after eq. (p^O]). We allow the exponent of the 
pomeron flux a to vary around a value (a = 1.2) which is smaller than the central value 
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(a = 1.3) of the ZEUS analysis [§] and closer to the most recent finding (a ~ 1.1 — 1.2) 
0. For comparison, the value used in Ref. is 1.20. For the effective cross section 
we take a^s = 10 mb. The x range shown in figure is the one allowed by the present 
experimental fits on the diffractive structure function. No decrease of shadowing occurs 
down to X ^ 10~^. 

In Fig. 2 we set a = 1.2 and we vary the coefficient C of the sea component of 



the diffractive structure function (10) within the errors of the ZEUS fit [^. We take 
(Teflf = 10 mb. Again, the shadowing curve is at most (when C is small) rather fiat 
towards x = 10~^ but does not exhibit any sensible decrease. 

Our main finding is therefore that the onset of shadowing saturation is at much 
lower X than argued in Ref. ||^ (a;o < lO"'*) and is outside the range of investigation of 
HERA experiments. Only a combination of unlikely circumstances (a and C very small, 
a < 1.1, C ~ 0, that is a pomeron fiux less singular than the Donnachie-Landshoff one 
and an almost vanishing sea component in the diffractive structure function) would 
produce a visible decrease of nuclear shadowing in the region above x = 10~^. 

In Fig. 3 we illustrate the dependence of our results on the effective cross section 
(Teff, for two nuclei (carbon and calcium). The shadowing curves with two choices of 
(Teflf (12 and 15 mb) and a = 1.2 in eq. ( |T0|) are shown. We see that the value of 
(Jeff affects the absolute normalization of nuclear shadowing but not its x-dependence. 
For light nuclei, such as carbon, not even the normalization depends on cXes- Thus the 
theoretical uncertainty related to the choice of cxeff does not spoil our conclusions about 
the X behavior of the shadowing curve. 

Considering the weak dependence of the shape of shadowing on the theoretical in- 
gredients of the present calculation (the effective cross section for multiple rescattering, 
the hadronic slopes and the longitudinal form factor), we can say that our prediction is 
dictated essentially by the experimental measurements of and therefore is a model 
independent result. Moreover, the qualitative behavior of shadowing that we found is 
stable against a large variation of the exponent of the pomeron flux and of the size of 
the sea diffractive structure function. Hence, any more precise determination of the 
diffractive structure function should not change our main flnding. 
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Finally, we just mention that results qualitatively similar to those reported above 
can be obtained by using F^^^^ evaluated in BFKL-type models, such as the one of 
Ref. |Tl|: only the absolute size of shadowing turns out to be larger. 

In conclusion, let us summarize our results. We studied the low-x behavior of 
nuclear shadowing at moderately large Q^. We carried out a model independent calcu- 
lation by using the experimental data on the diffractive structure function and the well 
established relation between diffraction and shadowing. No restrictive assumptions 
were made. Our conclusion is that, to the best of the present experimental knowledge, 
no decrease of shadowing occurs above x ~ 10"'*, that is in the kinematic region ac- 
cessible at HERA. Obviously, if the present experimental results on the different low-x 
behavior of F,^ and F^^^^ will be confirmed, a decrease of shadowing does take place, 
but only at very small x, beyond the HERA range. 



We thank M. Arneodo, E. Barberis, A. Solano and A. Staiano for many useful 
discussions on the ZEUS analyses. 
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Figure captions 



Fig. 1 The ratio F^{x)/F^{x) for calcium in the range 10-30 GeV^. The three 
curves correspond to three different values of the exponent of the pomeron flux 



in the ZEUS fit of the diffractive structure function eq. (|10D: a = 1.15 (dotted 
line), a = 1.20 (dashed line), a = 1.25 (solid line). The other parameters used 
are specified in the text. 

Fig. 2 The ratio F^{x)/F^ (x) for calcium in the range 10-30 GeV^. The solid 
curve is obtained using C = 0.57 in the ZEUS fit, see eq. ([T0|) . The dashed 
curves correspond to C = 0.57 ± 0.25 (errors added in quadrature). The other 
parameters are specified in the text. 

Fig. 3 The ratio F.f{x)/F,^{x) for carbon (upper pair of curves) and calcium (lower pair 
of curves) in the range 10-30 GeV^. The solid and dashed lines correspond to 
(Xcs = 10 mb and cTcfr = 15 mb, respectively. 
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